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ABSTRACT 

We compare the observed merger rate of galaxies over cosmic time and the frequency 
of collisional ring galaxies (CRGs), with analytic models and halo merger and collision 
rates from a large cosmological simulation. In the Lambda cold dark matter model 
(LCDM) model we find that the cosmic merger fraction does not evolve strongly 
between 0.2^z^2, implying that the observed decrease of the cosmic star formation 
rate since z~ 1 might not be tied to a disappearing population of major mergers. 
Haloes hosting massive galaxies undergo on average ~ 2 mergers from z^ 2 up to 
present day, reflecting the late assembly time for the massive systems and the related 
downsizing problem. The cosmic merger rate declines with redshift: at the present 
time it is a factor of 10 lower than at z~ 2, in reasonable agreement with the current 
available data. The rate of CRG formation derived from the interactions between halo 
progenitors up to z = 2 is found to be a good tracer of the cosmic merger rate. In the 
LCDM model the rate of CRGs as well as the merger rate do not scale as (1 + z) m , as 
suggested by previous models. Our predictions of cosmic merger and CRG rates may 
be applied to forthcoming surveys such as GOODS and zCOSMOS. 
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In a hierarchical Universe, galaxy mergers are thought to 
play an important role during structure formation, partic- 
ularly at higher redshifts. Mergers may be also relevant to 
the growth of massive early-type galaxies (e.g. Toomre 1977; 
Barnes & Hernquist 1996, Naab et al. 2006, Cox et al. 2006). 
Simulations of mergers involving gas-rich discs suggest that 
major mergers can trigger violent starbursts and transform 
discs into spheroidals. However the role of mergers in the 
galaxy assembly and star formation processes is still un- 
clear. The observed correlation between galaxy morphol- 
ogy and color indicates that the star formation history of a 
galaxy is closely tied to its morphology evolution. However, 
it has been recently recognized that there may be different 
timescales for the formation of stars in massive spheroidals. 
Thus, tracking the galaxy merger rate as a function of red- 
shift can constrain the contribution of mergers to the for- 
mation of stars in spheroids. 

Despite their importance, it has proved challenging to 
measure the rate of galaxy mergers and its evolution with 
cosmic time. Many theoretical and observational attempts 
have attempted to reconstruct the history of the galaxy in- 
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teraction rate (Toomre 1977; Zepf & Koo 1989; Carlberg 
1990a, 1990b; Carlberg, Pritchet & Infante 1994; Burkey 
et al. 1994; Yee & Ellington 1995; Neuschaefer et al. 1995; 
Woods, Fahlman, & Richer 1995; Patton et al. 1997; Le 
Fevre et al. 2000; Conselice, Bershady & Jangren 2000a; 
Conselice, Bershady & Gallagher 2000b; Bershady, Jangren 
& Conselice 2000; Conselice 2003; Conselice et al. 2003, 
2004; Cassata et al. 2005; Conselice 2006; Bridge et al. 2007; 
Jogee et al. 2007; 2008; see Conselice 2007 for a short re- 
view) . 

Various models (Toomre 1977; Carlberg 1990a, 1990b) 
suggest that the galaxy merger rate per unit volume n in- 
creases with the redshift z as: 



h oc (1 + z) r 



(1) 



The theoretical approach proposed by Carlberg (1990a, 
1990b), based on the Press— Schechter formalism (Press & 
Schechter 1974), predicts that the value of the coefficient m 
depends on the present-day matter density parameter Qm'- 



m ~ 4.51 x Q, ( m 2 . 



(2) 



A number of studies have used a sample of haloes or sub- 
haloes using N-body simulations to estimate merger rates as 
a function of redshift. Governato et al. (1999) studied major 
mergers of galaxy-sized haloes in open CDM universe. Got- 
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tloeber et al. (2001) studied the environmental dependence 
of major merger rates as a function of redshift in the con- 
cordance CDM models. More recently, the Millennium sim- 
ulation (Springel 2005) has been used to construct merger 
trees and to quantify the merger rates of haloes (Fakhouri 
& Ma 2008). The authors find that the average merger rate 
per halo depends weakly on the halo mass and that the halo 
merger rate evolves with z as (1 + z) m with m in the range 
2 to 2.3. Guo & White (2008) found similar results using 
the Millennium simulation galaxy catalog: the halo merger 
rates depend on redshift and only weakly on the halo mass. 
Murali et al. (2002) studied the relative contributions of 
merging and smooth accretion to the rate at which large 
galaxies gain mass using cosmological cubes simulated with 
N-body and SPH techniques. Mailer et al. (2006) estimated 
the merger rates of galaxies identified within subhaloes using 
cosmological hydrodynamic simulations. Additional studies 
used semi-analytic models to infer the merger fraction of 
haloes (e.g. Khochfar & Burkert 2001; Benson et al. 2005). 
Most of the models predict that the major merger rate of 
galaxy-sized dark matter haloes rises rapidly with redshift. 

The number of close pairs is often used as an obser- 
vational tracer of the galaxy merger rate. However obser- 
vational studies suggest of the number of close companions 
evolve less, over the cosmic time, than inferred from theo- 
retical studies of dark halo merging. Barrier et al. (2006) 
estimated the major merger rates of subhaloes using analyt- 
ical models plus numerical N-body simulations and studied 
the connection of the merger rates to the observed num- 
ber of close pairs. The authors use the halo-occupation- 
distribution (HOD) modelling of galaxy clustering to ex- 
plain the little evolution in redshift of the observed close- 
pair count. They find that the discrepancy is due to the 
additional processes occurring during merger of subhaloes 
in a common parent halo, which may be ignored in the halo 
merger rates. 

Due to the difficulty of directly observing the merger 
features, measurement of the redshift evolution of the frac- 
tion of galaxy pairs has traditionally been taken, which can 
be parametrized as oc (1 + z) . Most studies (Zepf & Koo 
1989; Carlberg et al. 1994; Burkey et al. 1994; Yee & Elling- 
ton 1995; Woods et al. 1995; Patton et al. 1997; Kampczyk 
2007) derive a value of k ~ 2 — 4, while Neuschaefer et al. 
(1995) find k ~ 0, due to a different estimate of the non- 
physical galaxy pairs. The main difficulty of this method is 
that the conversion from k to m is unclear. It is also diffi- 
cult to disentangle projection pairs from true physical pairs 
when using photometric redshifts alone. Additionally, not all 
galaxies in a physical pair will merge, as the galaxies may 
be unbound. 

From an observational point of view, it is difficult to es- 
timate whether h depends on the redshift and to determine 
the correct value of m. A direct measurement of merger fea- 
tures in distant galaxies is difficult, as tidal tails and distor- 
tions generally have low surface brightness (see Mihos 1995; 
Hibbard & Vacca 1997). 

A powerful method for measuring the galaxy merger 
rate is to count the incidence of strongly disturbed galaxies 
(with strong asymmetries, double nuclei or prominent tidal 
tails), the so-called CAS method (Conselice 2003). There are 
several methods to quantify the frequency of strongly dis- 
torted galaxies: visual classification, quantitative measures 



of asymmetries such as the CAS system (Conselice 2003) 
and the Gini-M20 system (Lotz et al 2006). 

All asymmetries suffer from surface brightness (SB) 
dimming, but the outer low SB features suffer more strongly 
from it. CAS misses the latter features, but visual classifica- 
tion captures many of these. Simulations (Conselice 2006) as 
well as empirical studies (Jogee et al. 2007, 2008) show that 
visual classifications capture a larger fraction of strongly dis- 
torted galaxies than the CAS merger criteria, as the eye is 
sensitive to asymmetries over a larger dynamic range. 

The CAS method, applied to the Hubble Deep Field 
(HDF), provides an estimate of m ranging from 4 to 6 (Con- 
selice et al. 2003; but m ~ 2 — 4 in the re-analysis of these 
data by Conselice 2006). Based on the results of the CAS 
analysis, Conselice (2006) suggests that equation Q is inac- 
curate for some galaxy types (especially for small galaxies) 
and for high redshifts (z £ 1 — 2). The CAS method provides 
a more straight-forward estimate of merger rate than the 
incidence of close pairs of galaxies. However, positioning a 
galaxy in the CAS plane is sometimes problematic. Further- 
more, the connection between high asymmetry/lumpiness's 
and merger history implies a number of assumptions. 

For these reasons, Lavery et al. (2004; hereafter L04) 
proposed to use ring galaxies as a more direct tracer of 
galaxy mergers. In fact, a high fraction of ring galaxies 
[« 60 per cent, Few & Madore (1986)], called 'P-type ring 
galaxies', are thought to have collisional origin. Recent N- 
body simulations (Mapelli et al. 2008a, 2008b, and refer- 
ences therein) show that the ring phase is quite short-lived: 
it lasts only for ^ 500 Myr after the galaxy collision. More- 
over, ring galaxies are easier to identify than other interac- 
tion signatures (e.g. tidal tails; see Mihos 1995; Hibbard & 
Vacca 1997). Thus, the number of collisional ring galaxies 
(hereafter CRGs) may be a straight-forward tracer of the 
galaxy interaction rate. 

Unfortunately, most ring galaxies with measured dis- 
tances are relatively nearby [z 0.1; see e.g. the sample of 
68 ring galaxies in Few & Madore (1986)]. 

L04 analyze 162 Wide Field Photo Camera 2 (WFPC2) 
fields, obtained from the Hubble Space Telescope (HST) 
Archives, in order to identify distant CRGs. They find 25 
CRGs in their images. From this sample, L04 derive a value 
of the merger rate m ~ 5. However, their estimate is affected 
by large uncertainties, as they have redshift measurements 
only for 6 of their 25 CRGs. For the remaining 19 galax- 
ies, they derive an 'estimated redshift', by assuming that 
CRGs have similar visual magnitude. Recently, Elmegreen 
& Elmegreen (2006, hereafter E06) analyzed other 24 CRGs 
in the GEMS and GOODS fields. For these galaxies redshift 
measurements are available. 

In this paper we present an attempt to derive the merger 
and the CRG formation rate from cosmological simulations. 

The simulations can estimate the rate of minor and ma- 
jor mergers in progenitors of the present-day galaxy haloes 
and establish whether the merger rates and CRG formation 
rates are related. We calculate the evolution of the merger 
and CRG formation rate up to redshift 2. We use a cosmo- 
logical cube a factor of three larger and higher numerical 
resolution than adopted in the SPH simulations of Mailer et 
al. (2006) and our numerical resolution is higher of almost 
a factor of ten than the Millennium run used in Fakhouri & 
Ma 2008). Finally this study compares results from numeri- 
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cal simulations with the available data and gives predictions 
for future observations. In Section 2 we present details of 
the numerical simulation and analysis procedure. Section 3 
discusses our main results, while Section 4 summarizes our 
conclusions and implications. 



2 NUMERICAL METHODS 
2.1 Simulations 

We analyze a cosmological iV-body simulation of the 
Lambda cold dark matter (LCDM) cosmogony, with cosmo- 
logical parameters chosen to match the 3-yr Wilkinson Mi- 
crowave Anisotropy Probe (WMAP3) constraints (Spergel 
et al. 2007). These are characterized by the present-day mat- 
ter density parameter, Qm = 0.238, a cosmological constant 
contribution, £7a=0.762, and a Hubble parameter h = 0.73 
{Ho = 100 h km s" 1 Mpc" 1 ). The mass perturbation spec- 
trum has a spectral index of n = 0.951, and is normalized 
by the linear rms fluctuation on 8 Mpc/h radius spheres, 
as = 0.75. 

We follow the evolution of 600 3 particles of mass rridm = 
8.67 x 10 7 h' 1 M e in a box of 90 Mpc (or 65.7 Mpc/h) (co- 
moving) on a side, by using the TV-body code PKDGRAV 
(Stadel 2001). Gravitational interactions between pairs of 
particles are softened with a fixed comoving softening length 
of 1.16 kpc. 



2.2 Halo Identification at z=0 

Non-linear structures at z = are identified using the clas- 
sic friends-of-friends (FOF) algorithm with a linking length 
equal to 0.2 times the mean comoving interparticle separa- 
tion. For each FOF halo we identify the most bound particle 
and adopt its position as the halo centre. Using this centre, 
we compute the 'virial radius' of each halo, r v i r , defined as 
the radius of a sphere of overdensity A(z = 0) = 94 (rel- 
ative to the critical density for closurejj. Quantities mea- 
sured within r v i r will be referred to as 'virial', for short. 
We select for our analysis all haloes with masses in the 
range M vir = 5 x 10 12 to 10 14 h' 1 M©. The resulting haloes 
have Avir between 56,000 and 1,100,000 particles within 
the virial radius. 



2.3 Merger tree Construction 

For each of the haloes in our z — sample we have con- 
structed a merger tree over the period < z < 2, based on 
FOF haloes. We use the FOF merger tree to define and quan- 
tify the accretion and merger rates of haloes. Other authors 
have previously used this technique or used substructure di- 
rectly. We prefer to avoid the use of substructure since at 
a fixed mass resolution, smaller haloes are less resolved and 
have less substructure than larger haloes due to the classic 
overmerging problem (Moore et al 1996). 



1 The virial overdensity in a flat universe may be computed using 
the fitting formula proposed by Bryan & Norman (1998): A(z) = 



187T 2 +82/(z) -39/(z) 2 ; with f(z) 



Oo(i+zr 
Q (i+z)3+n A 



- 1 



We consider a FOF halo identified at z > to be 'pro- 
genitors' of a z = system if at least 50 per cent of its par- 
ticles are found within the latter. Using this definition we 
can identify, at all times, the list of progenitors of a given 
z = halo and track their properties through time. In the 
tree to identify halo mergers, we denote a halo as a major 
merger remnant if at some time during <z< 2 its major 
progenitor was classified as a single group in one output but 
two separate groups with a mass ratio ^ 4 : 1 in the pre- 
ceding output (see D'Onghia & Burkert 2004). As already 
noted in previous works (e.g. Gottloeber et al. 2001, Berrier 
et al. 2006, Fakhouri & Ma 2008), the merger tree can result 
in fragmentation events, in which particles of a progenitor 
halo end up in two distinct halo descendants. This spurious 
fragmentation is an artifact of the FOF halo identification 
scheme. During the initial merge phase the halo finder asso- 
ciates and dissociates particles inside and outside the bound 
region. We noted that the fragmentation event of a halo 
progenitor lasts for less than 0.5 Gyr. After this time usu- 
ally the progenitor merges again. If this fragmentation is not 
properly treated, the risk is to count twice the same merger 
event for the current halo. We avoid this risk by assuming 
that the merger happens the first time the progenitor halo 
is considered part of the descendant. Our method to handle 
fragmentation events is similar to the 'stitching' method as- 
sumed from Fakhouri & Ma (2008). Each FOF halo at z > 
in the tree catalog is assigned a mass counting the number of 
particles associated with the FOF group, rather than a mass 
defined with the overdensity criterion. Halo progenitors of 
interest contain at least 250 particles. 

An additional problem with the merger tree is that the 
mass of a descendant halo is not exactly the sum of all the 
progenitors, but it may count inside the halo a diffuse mass 
component which is not resolved in subhaloes due to lim- 
ited numerical resolution. The diffuse mass component may 
dominate the merger events at high redshift, when it is more 
difficult to resolve progenitors, owed the limited numerical 
resolution. However our current numerical resolution is al- 
most a factor of ten higher than the Millennium run used 
from Fakhouri & Ma (2008) and of the SPH simulation of 
Mailer et al. (2006) and allows us to resolve progenitors and 
relative mergers and guarantees that we do not miss any 
relevant major merger event and interaction of haloes pro- 
genitors back in time. 

Previous studies built the halo merger trees by connect- 
ing subhaloes instead of FOF haloes across the snapshot 
outputs (Gottloeber et al. 2001, Mailer et al. 2006, Berrier 
et al. 2006). In these papers the procedure to define de- 
scendants and progenitors is similar to that one assumed for 
FOF trees: a subhalo at a given redshift is assumed to be de- 
scendant of a progenitor subhalo defined at higher redshift 
if it contains a fixed percentage of particles. In particular 
Berrier at al. (2006) uses a hybrid N-body simulation plus 
analytic substructure model to predict the number of pairs, 
which is the quantity often used to infer the observed galaxy 
merger rate. Since the observed number of close companions 
rises with redshift slower than the halo major merger rates 
predicted from previous simulations, they assume a halo- 
occupation model, in order to match the observed number 
of close companions inferred from current data. Indeed, our 
algorithm is based on a FOF halo tree and does not include 
subhaloes. However, a merger tree constructed on subhaloes 
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inside a parent halo is mainly required to analyze galaxy- 
cluster sized haloes and massive galaxy-group sized haloes. 
The cosmological cube considered here contains only one 
cluster-sized halo. The galaxy-sized haloes and galaxy-group 
sized haloes defined in our sample are located in low density 
environments and are compared to galaxies which are not 
located in clusters, but mostly in the field. 

It worth noting that when merger trees are based on 
subhaloes it is difficult to estimate the mass of a subhalo 
located within a larger halo. Tidal stripping occurring when 
a subhalo enters into larger halo can reduce significantly 
the mass of the subhalo before it approaches the pericentric 
distance, with implications for the inferred mass ratio of the 
merging events. Since our merger trees are based on FOF 
haloes and do not consider subhaloes, our analysis does not 
suffer of this problem. 

2.4 Ring Galaxy Formation Criteria 

Since our simulations are based on dark matter only, we 
cannot directly trace the formation of CRGs. However, we 
can estimate the rate of CRG formation from our models 
adopting the following criteria. At each redshift we select 
a sample of progenitor haloes in the range of mass of a few 
10 11 Mq to 5 x 10 12 Mq (corresponding to the typical masses 
of observed CRGs) and we consider the encounters between 
them. We consider a list of 8 progenitors ordered by de- 
creasing mass and we examine the collisions between the 
most massive and the second most massive progenitor and 
all the possible combinations amongst progenitors along the 
list. 

We assume that a CRG is formed whenever two haloes 
undergo an encounter in which: 

i) the mass ratio between the bullet (hereafter 'intruder' 
halo) and the most massive progenitor (hereafter 'target' 
halo) is > 1 : 10; 

ii) the pericentric distance p is less than ~ 15 per cent of 
the expected disc radius of the target galaxy (i.e. p ^ 4 kpc, 
Lynds & Toomre 1976). 

The first constraint comes both from observations of nearby 
CRGs whose intruder is known, and from numerical simula- 
tions showing that the ring is hard to form when the mass 
ratio between the intruder and the target galaxy is ^ 1 : 10 
(Hernquist & Weil 1993; Mihos & Hernquist 1994; Horellou 
& Combes 2001; Mapelli et al. 2008a, 2008b). 

The second condition assumes results from numerical 
simulations showing that circular rings form only when the 
impact parameter is small. The larger the impact parame- 
ters, the more asymmetric is the resulting ring (Hernquist 
& Weil 1993; Mihos & Hernquist 1994; Horellou & Combes 
2001; Mapelli et al. 2008a, 2008b). 

Finally, we do not put any limitation on the inclination 
angle 9 (between the disc axis of the target and the velocity 
of the intruder). Lynds & Toomre (1976) indicated that rela- 
tively symmetric ring galaxies can form at least for 9 ^ 45° . 
Recent simulations (Ghosh & Mapelli 2008) show that reg- 
ular (although warped) rings form also for 9 > 60°. So, we 
can reasonably assume that ring galaxies can form also for 
high values of 9. Thus, our estimate of the CRG formation 
rate represents an upper limit and should be rescaled by an 
unknown factor 1 — 003(6*) ^ 0.5. 



3 RESULTS 

3.1 Merger Fraction Evolution up to Redshift 2 

To predict the merger fraction as a function of redshift using 
observations of galaxy pairs is a difficult task, due to the 
complexity of establishing the time over which the merger 
occurred. Estimates of merger fractions using galaxy pairs 
come from Patton et al. (2000, 2002), Le Fevre et al. (2000) 
and Lin et al. (2004) for z < 1. However, owing to projection 
effects, some pairs might not be a physically bound system 
and merging may last for a long time or may never occur. 

On the other hand, in simulations we assume that a 
merger occurs when the smaller halo enters into the virial 
radius of the host halo, whereas the distance between the 
observed pairs is generally smaller than the virial radius of 
the progenitor. In fact, the baryonic component of the two 
galaxies is within ~ 20 per cent of the virial radius. 

A comparison of the merger rates extracted from the 
simulations for dark haloes to galaxy merger rates has some 
limitations. A direct comparison implies a conversion be- 
tween halo mass and galaxy mass/light. Many models usu- 
ally adopted to this conversion, e.g., see van den Bosch et al. 
(2007) and reference therein, have a strong dependence upon 
halo (or galaxy) mass. These models imply that a merger be- 
tween a dark matter halo and another halo of one-tenth its 
mass may not be equivalent to a merger between a galaxy 
and another one of one-tenth its mass. Over the mass range 
10 n A'/ Q to 10 12 M Q , the variation in M/L ratio with halo 
mass is approximately a constant in the HOD models, so 
this approximation is good. For larger mass haloes we should 
compress our merger mass scale by 30%, but this depepence 
is model based and would introduce an error smaller than 
the observations currently suffer from. Ideally we would use 
hydrodynamical simulations that accurately resolve the star- 
formation and luminosity evolution of galaxies. This is a few 
years away before becoming achievable computationally. 

The timescale over which the merger between a pair 
of galaxies is deemed visible is estimated as a fraction of 
the crossing time t CI of the infalling halo into the host halo. 
The crossing time for realization of 1:1 merger simulations 
is ~ 1 Gyr. Simulations of 1:1-1:2-1:3 mergers show that 
the timescale over which merger features are visible is ~ 0.5 
Gyr for systems with different orbital inclinations (Conselice 
2006). 

We therefore define the halo merger fraction / mer g as 



,( f> Af) = 

= v ; iVtot At' 



(3) 



number of haloes undergoing a (minor or major) merger at a 
given time interval At and for a given mass M, respectively. 
t CT /At is the fraction of the time interval over which the 
merger signature is visible. In the time interval between z=0 
and z=0.3 the time is longer than the merger time so we did 
not account for the time delay t CT / At. 

Table 1 lists the number of major and minor mergers 
and the total number of haloes derived in simulations. In 
Figure [T] major merger fractions (the dotted line) and total 
(i.e. minor+major) merger fractions (solid line) derived from 
simulations are plotted as a function of redshift when assum- 
ing t cr = 1 Gyr. Our models predict a merger fraction that 
does not evolve significantly with redshift between 0.2^z^2. 
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Table 1. Minor and Major Mergers. 



z 


Mergers 


Major Mergers 


b 

Mergers 
[M~ 10 12 M Q ] 


Major Mergers 
[M~ 10 12 M Q ] 


IN tot ildlOCS 


0.0 


571 


11 








0.1 


837 


19 


3 


3 


756 


0.3 


649 


30 


11 


7 


1565 


0.5 


655 


17 


30 


8 


2232 


0.7 


706 


43 


65 


26 


2844 


0.9 


496 


4 


71 


3 


3436 


1.1 


424 


21 


91 


12 


3852 


1.3 


326 


6 


87 


2 


4240 


1.4 


254 


10 


91 


8 


4493 


1.6 


267 


15 


88 


7 


4669 


1.7 


302 


14 


82 


7 


4813 


1.9 


210 


14 


69 


6 


5020 



a Total Mergers (minor+maj or) and major mergers from simulations for haloes at present day in the range of mass between 5 X 10 12 
and 10 14 Mq. 

Total Mergers(minor+major) and major mergers from simulations for progenitors with fixed mass between 9 X 10 11 and 3 X 10 12 Mq 
at any redshift. 
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Figure 1. Total merger fraction (solid line) and major merger 
fraction (dotted line) of dark haloes with mass between 5 X 10 12 
and 10 14 Mq at present day, as a function of redshift, derived from 
cosmological simulations. The long-dashed line shows the major 
merger fraction of Milky- Way sized haloes at any time. Symbols 
represent empirical results on the merger fraction or fraction of 
strongly distorted galaxies identified via different methods, rang- 
ing from visual classification to automated methods: Jogee et al. 
(2007, 2008, filled circles), Le Fevre (2000, filled squares), Con- 
selice (2003, open squares) and Lotz et al. (2006, stars). Note 
that a stellar mass cut off was applied in Jogee data points 
{M, > 2.5 X 10 10 Af Q ) and in Conselice 2007 (M* > 1 X 10 10 M Q ). 



Wolf et al. (2005) and Bell et al. (2005) atz~ 0.7. Table 1 
also lists the number of total mergers (minor+major) and 
major mergers for Milky Way sized progenitors identified at 
any time, with typical mass of ~ 10 12 Mq. 

In Figure [T] we compare the predictions of our models 
with data of the merger fraction, or fraction of strongly dis- 
torted galaxies, identified via different methods and based on 
different surveys : Jogee et al. (2007, 2008, filled circles), Le 
Fevre (2000, filled squares), Conselice (2003, open squares) 
and Lotz et al (2006, stars). The results of Jogee et al. (2007, 
2008) refer to the the fraction of strongly distorted interact- 
ing/merging massive (with stellar mass M* > 2.5xlO lo M ) 
galaxies over z ~ 0.24 to 0.80, identified from the GEMS sur- 
vey (Rix et al. 2004) using both the CAS system and an in- 
dependent visual classification system, specifically designed 
to separate interacting galaxies with externally-triggered 
asymmetries from non-interacting galaxies with small-scale 
internally-triggered asymmetries . 

The presence of wiggles in the theoretical estimates is 
due to the uncertainty of using the merger tree to identi- 
fying mergers between progenitors at each time. The major 
merger fraction of Milky- Way sized haloes identified at any 
time is displayed with the long dashed line. We note that the 
predicted fraction of major mergers and minor mergers is al- 
most constant from z=2 up to present day for Milky- Way 
sized haloes identified at any redshift and for all the haloes 
of our sample, regardless the mass. Furthermore, the agree- 
ment between merger fractions predicted from cosmological 
simulations (based on the halo merger history) and the ob- 
served merger fractions based on galaxy CAS morphologies 
is encouraging for the LCDM model. 



This implies that the observed decrease of the cosmic star 
formation rate since z~l (e.g. Lilly et al. 1996; Madau et 
al. 1996) is not tied to a disappearing population of major 
mergers, and seems to be in agreement with recent new re- 
sults by Jogee et al. (2007; 2008) for z ~ 0.2 - 0.8, and by 



3.2 Merger Rates 

The merger fraction is related to the merger rate per unit 
volume 7?, defined within a time interval and mass range, 
by the following expression. 

R(t,M) = i merg T~ 1 n m , (4) 
where r m is the timescale for a merger to occur and 
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Figure 2. Total merger rate (major+minor) (solid line) and ma- 
jor merger rate only (dotted line) of dark matter haloes, per unit 
volume (comoving), as a function of redshift, derived from cos- 
mological simulations. The rate per unit volume of Milky Way 
sized haloes identified at any time is plotted for comparison (long- 
dashed line, blue on the web). Filled circles: rate of mergers and 
interactions, with mass ratios in the range 1:1 to 1:10 , derived 
by Jogee et al. (2007, 2008) from the GEMS data. Open squares: 
merger rate from the HDF from Conselice (2003). 



is defined as ~ t cr and n m is the physical density of the 
haloes undergoing a merger (minor or major) within a given 
mass range and at a given time. The physical density of 
merging haloes is derived by dividing the total number of 
haloes which are merging N meis , listed in Table 1, by the 
physical volume occupied by all the considered haloes. 

Figure [2] plots the major merger rate (dotted line) and 
the total merger rate (minor+major) (solid line) per unit 
volume, as a function of redshift, as derived from cosmolog- 
ical simulations. For comparison, Figure [2] also shows the 
merger rate inferred from the GEMS sample (Jogee, private 
communication; marked with filled circles) and the analy- 
sis of Conselice et al. (2007) based on the Hubble Ultra Deep 
Fteld (HUDF, open squares). For completeness, the plot dis- 
plays the major merger rate per unit volume inferred in cos- 
mological models for progenitors in the range of mass of the 
Milky way (~ 1O 12 M0), identified at any redshift (long- 
dashed line). A large fraction of the progenitors are in this 
range of mass at higher redshift, explaining why these sys- 
tems show a similar trend in the merger rate as the major 
mergers of all the haloes. The merger rate of Milky Way 
sized haloes (long-dashed line) decouples from the total ma- 
jor merger rate (dotted line) only at z 0.5, where the for- 
mer drops, while the latter decreases more gently. This is 
due to the fact that mergers occurring at late time mainly 
involve the assembly of larger mass systems. 

Note that the decrease of the volume-averaged merger 
rate at late times is a result of the decrease of the progen- 
itor number density at lower redshift. Our predictions are 
in good agreement with the observations, despite the large 



uncertainties in the available data. It worth noting that our 
estimates based on dark matter haloes are in agreement with 
the estimates inferred from Mailer et al. (2006) who used N- 
body simulations plus hydrodynamics with a merger tree 
based on subhaloes to better trace the close pairs of galax- 
ies. It is encouraging to note that estimates based on dark 
haloes and SPH simulations do agree without invoking any 
additional treatments of baryons in galaxies with the halo- 
occupation models and semi-analytic models. 

Summing over all mergers for massive haloes (> 5 x 
1O 12 M0), we find that the average number of mergers a 
halo experienced since z ~ 2 is N m ~ 2 in agreement with 
estimates of Conselice et al. (2007). 



3.3 Ring Galaxy Formation Rate 

First we check the correspondence between the CRG forma- 
tion rate and the halo merger rate within our cosmological 
simulation. For CRGs we consider only halo progenitors with 
mass between a few 10 11 and 5 x 10 12 Mq at any time, which 
correspond to the observed masses of ring galaxies. Figure [3] 
shows the CRG formation rate (lightly hatched histogram, 
red on the web), compared with the total (minor+major) 
merger rate of all the haloes with present day masses be- 
tween 5 x 10 12 Mq and 10 14 Mq (open histogram, blue on 
the web). Both the CRG formation rate and the merger rate 
increase at lower redshifts. The null hypothesis probability 
that ring galaxies and mergers are drawn from the same 
distribution from redshift 0.2 up to z=2 is ~ 0.11, corre- 
sponding to a non-reduced \ 2 = 17.1 (for 11 data points, 
parameters, and assuming Poissonian errors). Thus, we 
conclude that the CRG formation rate can be considered a 
tracer of the merger rate. Since our results rely upon the dis- 
tribution of merging orbits we checked the consistency with 
previous works. We computed the fraction of merging pro- 
genitors with different pericentric distances using the sample 
we adopted for the CRG calculations. By tracking accretion 
events we predict the pericentric distance of the merging 
haloes. This distribution does not disagree with that shown 
in Khochfar & Burkert (2006) and is similar to that found in 
simulations that resolve substructures in CDM haloes (e.g. 
Ghigna et al 1998). 

Finally, we remark that our simulations follow the 
merger history and dynamics of dark haloes and neglect 
processes involving the baryonic physics. In particular, the 
formation of galaxy discs, which is crucial for the formation 
of CRGs, cannot be modeled here. Also the time evolution 
for the late and early type number fraction is not accounted. 



3.4 Comparison with Observations 

The available data of moderately high redshift CRGs consist 
of a sample of 25 ring galaxies (0.2 ^ z ^ 1), observed with 
HST (L04), and a more recent sample of 24 ring galaxies 
(0.07 < z «S 1.5) found in the GEMS and GOODS fields 
(E06). Table 2 reports the number of CRGs observed in 
these two samples, as well as the number derived from our 
simulations. 

Unfortunately, a direct redshift estimate has been pro- 
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Figure 3. CRG formation rate (i.e. the number of newly formed 
CRGs in the entire simulation per redshift interval, lightly 
hatched histogram, red on the web) as a function of redshift, com- 
pared with the merger rate obtained considering all the haloes in 
the range of mass between 5 X 10 12 Mq and 10 14 Mq at z = 
(open histogram, blue on the web; see Table 1 and Fig. 2). The 
histograms are normalized to total number of CRGs and mergers, 
respectively 



Table 2. Number of simulated and observed CRGs. 



Figure 4. CRG formation rate per unit volume (comoving). Open 
triangles (case L04old, green in the online version): the sample of 
25 observed CRGs in L04 with z es t (the same as in figure 4 of 
L04). Filled squares (case L04z o (, s , blue in the online version): 
the 8 observed CRGs with z obs (CRG number 2, 3, 5, 7, 9, 10, 12 
and 20 of L04). Filled triangles (case L04new, black in the online 
version): the 8 observed CRGs with z oos and the remaining 17 
observed CRGs with Zest- Open circles (magenta in the online 
version): CRG rate per unit volume derived from the E06 data. 
Long dashed line (red in the online version): simulated CRGs. 
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a CRGs from our simulation (considering haloes with mass 
between ~ 10 11 and 5 X 10 12 Mq). The total number of CRGs 
in the simulation, from z = 0.1 to z = 2 (from z = 0.2 to z = 1) 
is 57 (23). 

b The redshift is observed for 8 galaxies in the L04 sample (see 
the Text for details). An estimated redshift (L04) has been used 
for the 17 CRGs without redshift measurement. Note that the 
area of the survey is different from the one of our simulation 
(see text for details). 

vided only for 8 of the CRGs of the L04 samplfl LQ4 pro- 



2 Redshift estimates are provided for 6 of the 25 galaxies of the 
L04 sample, in particular the estimates reported in Table 1 of 
L04 refer to CRGs identified with number 2; 3; 5; 7; 10; 20. 
Two more galaxies of that sample (CRGs labeled as 9 and 12) 



posed to estimate the redshift of CRGs by adopting a 'stan- 
dard' absolute V magnitude for CRGs (Appleton & Marston 
1997). 

Note that only the estimated redshifts (z est ) are shown 
in figure 4 of L04, even for those CRGs that have an observed 
redshift (z bs). Also, the difference between z es t and z a bs is 
?L 10 per cent, i.e. half of the bin width, for 4 of the 8 
CRGs with known redshift. Yet, for 3 of these 4 galaxies 
the observed redshift z bs is significantly smaller than the 
estimated value z eat ■ 

This discrepancy has important consequences in the es- 
timate of the CRG formation rate. In Figure |4]the CRG for- 
mation rate per unit volum«[f| of the CRG sample observed 
by L04 is shown as a function of redshift. For the open 
triangles (hereafter case 'L04old') of Figure [4] we assumed 
the estimated redshift z est reported in L04 for all the 25 
CRGs (including the 8 CRGs with measured redshift z b 3 ). 
Note that the case L04old is the same as figure 4 of L04. 
The formation rate of the 8 CRGs with observed redshift 
z bs is represented by the filled squares of Figure [4] (here- 
after 'L04z o f, s '). Note that the ring galaxy formation rate of 
the sample of L04 peaks at z = 0.7 when z es t is assumed 



are members of galaxy clusters CL 0303+1706 (z = 0.6564) and 
CL 1601+4253 (z = 0.5382), respectively (see e.g. Dressier & 
Gunn (1992) for the redshift determination). 

3 The total solid angle of the observations in L04 and in E06 was 
~ 7.31 X 10 -5 sr and ~ 4.87 X 10 — 5 sr, respectively, independent 
of redshift. 
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(L04old), whereas the peak shifts to z — 0.5 when z t s is 
assumed instead (L04z O 6 S ). 

For completeness, we plot the rate of formation of 25 
CRGs of L04 when z bs is assumed for 8 galaxies which have 
an observed redshift and z est assumed for the remaining 17 
CRGs (filled triangles, hereafter 'L04new') . The distribution 
peaks at z ~ 0.7. The open circles in Figure represent 
the CRG formation rate derived using the data reported by 
E06. In this case all the redshifts are measured (either with 
spectroscopy or photometry). 

We additionally plot in Figure [4] the CRG formation 
rate per unit volume derived from cosmological simulations, 
as a function of redshift (dashed line). Note that the CRG 
formation rate is derived from Table 2 and is the same 
as shown in Figure [3] but with a different normalization. 
The simulated CRG formation rate per unit volume ap- 
proximately matches the data points in the redshift range 
0.2 sC z < 0.8. 

We also note that in our simulations the CRG formation 
rate per unit volume substantially increases with redshift. 
As already noted for the merger rate, this is mainly due to 
a substantial increase of the number density of progenitors 
present in merger tree at higher redshift. 

At z ~ 1.1 — 1.5 our simulated CRG formation rate 
is a factor of ~ 6 — 22 higher than the one derived from 
E06 data. This might be due to several factors. First, we 
assume that all progenitors with mass 10 11 — 5x 10 12 Mq and 
experiencing interactions with small impact parameter are 
disc galaxies, but the morphology of the progenitor cannot 
be inferred from our dark matter only simulations. Second, 
the E06 sample might be considered incomplete at z > 1.1. 

Thus, new redshift measurements of moderately high- 
redshift CRGs will be extremely useful. The future survey 
zCOSMOS (Lilly et al. 2007) will acquire spectra and red- 
shifts of approximately 10,000 galaxies (0 < z < 3) in the 
COSMOS survey field and will derive for 10,000 galaxies be- 
tween <z< 1 the mass, the morphology and the size. This 
will provide an excellent sample to compare with the rate of 
CRGs predicted from our models. 

3.5 Evolution of the Number of CRGs with 
Redshift 

Various theoretical models suggest that the galaxy merger 
rate per unit volume scales with the redshift as h cx (l + z) m . 
The exact value of m depends on the details of the adopted 
formalism, as well as on the cosmological parameters. Ap- 
proaches based on the Press-Schechter formalism (Carlberg 
1990a, 1990b) predict m ~ 2.5, assuming Qm = 0.238 from 
the WMAP3 constraints. 

We showed in Figure [3] that the CRG formation rate 
the is a good tracer of the merger rate. Hence, one might 
assume that also the CRG rate scales with (1 + z) m (L04). 
Thus, we derive the value of m from the simulated CRGs. 
In particular, the expected number of CRGs (JVcr,g(zi, 22)), 
which form between Z\ and zi, in a given volume, under the 
assumption that the density of CRGs scales as in equation]]] 
may be expressed by the following formula (see section 3.2 
of L04). 

Nckg(z!, z 2 ) = ncRGfi (;#-) X 



Table 3. Number of expected CRGs in our simulation, as a func- 
tion of m, in the redshift range zi, zi = (0.2, 1.0) (central col- 
umn) and 21, Z2 = (0.2, 2.0) (right column). 
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where c is the light speed, ticrg,o ~ 5.4 x 
10~ 6 h 3 Mpc~ 3 is the current density of CRGs (Few & 
Madore 1986), £{z) = [(1 + z) 3 Q M + 0\] 1/2 and Afi(z) 
is the considered solid angle at a given redshift z. 

In the simulations Afl(z) is the total solid angle, at 
a given z, occupied by our halo sample. Att(z) has been 
calculated as the sum of the physical sizes of each halo in 
the simulation divided by the comoving distance at a given 
z, i.e. 
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where a is the cosmological factor of expansion. The equa- 
tion ([6| is normalized to 2 = 0.3. 

Solving equation (JS| for the values of Afl(z) derived 
in equation (|6} gives the expected number of CRGs which 
form in our simulations between zi and Z2, provided that 
the CRG density scales with (1 + z) m . 

We calculated equation J5} for different values of m and 
for (zi, Z2) = 0.2, 1.0 and (21, 22) = 0.2, 2.0. Results are 
reported in Table 3. Our cosmological simulations show that 
23 CRG form between Z\ = 0.2 and 22 = 1.0, leading to the 
value of m ~ 3.7, which is lower than the one obtained by 
L04 (m ~ 5). However, if we consider the simulated CRGs 
forming up to 22 = 2, m is slightly lower: 46 CRGs are 
expected to form if m ~ 3.3. 

If we adopt the model by Carlberg (1990a, 1990b), a 
value of m = 3 — 4 is in disagreement with the cosmologi- 
cal parameters measured by WMAP3. In fact, equation ((2| 
implies Qm ~ 0.38 — 0.76 for m = 3 — 4, at odds with the 
estimate of WMAP3. 

Viceversa if Qm = 0.238 is assumed, in agreement with 
WMAP3 data, a value of m ~ 2.5 is derived from equa- 
tion @. For m ~ 2.5 the predicted number of newly formed 
CRGs in simulations up to 2 ~ 1 is ~ 13 (see Table 3), too 
low with respect to the observations of L04. 

Furthermore, in the bottom panel of Figure [S] we plot 
the evolution in redshift of the number of CRGs obtained 
in cosmological simulations (hatched histogram) compared 
with the number of CRGs predicted from equation © (open 
histogram with dashed line). We note the different evolution 
of the two distributions with redshift: the former decreases, 
the latter increases with redshift. The null hypothesis proba- 
bility that the two histograms are drawn from the same dis- 
tribution is 0.19, corresponding to a non-reduced x 2 = 12.4. 
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Figure 5. Bottom panel: number of newly formed CRGs as a 
function of redshift. Hatched histogram (red in the online ver- 
sion): simulated CRGs. Open histogram with dashed line: model 
expressed by equation JSJl for m = 3.3 (see the Text and L04). 
Open histogram with dotted line (blue on the web): model ex- 
pressed by equation (0 for 7 = 50, m = 2.4 and j3 = —2 [see 
the Text and Conselice (2006)]. Open histogram with dot-dashed 
line (green on the web): flat distribution. Top panel: total merger 
rate (hatched histogram, red on the web). Open histogram with 
dashed line: equation for m = 3.3. Open histogram with dot- 
ted line (blue on the web): equation J7J for 7 = 5000, m = 2.4 
and (3 = -2. 



(with 10 data points, 1 parameter, and assuming Poissonian 
errors) . 

We also compare the number of CRGs obtained in our 
cosmological simulations with the fitting function proposed 
by Conselice (2006) [see equation (5) of Conselice (2006)]. 
The fitting formula by Conselice (2006), adapted for our 
simulated CRGs, can be written as 

Ncrg(z) = -y(l + z) m exp[/3(l + z)], (7) 
where 7=50, m=2.4 and /3=-2.0 best fit the simulated 
CRG distribution (between z — 0.2 and z — 2.0), reported 
in Figure [5] 

In this case, the non-reduced x 2 is 7.5, leading to a null 
hypothesis probability equal to 0.38 (considering 10 data 
points, 3 parameters, and assuming Poissonian errors). 

However, the simulated CRG rate can be fit also by a 
flat distribution with TVcrg = 5.0 (non-reduced \ 2 = 12.4 
and null hypothesis probability equal to 0.19, for 10 data 
points, 1 parameter and Poissonian errors). 

In Section 3.3 we stressed that the CRG formation rate 
is a good tracer of the merger rate. For comparison, the top 
panel of Figure [5] shows the evolution with redshift of the 
total merger rate (i.e. the merger rate of haloes with mass 
between 5 x 10 12 and 10 14 M at z = 0). 

In this Figure, the total merger rate is compared with 
the Carlberg model (dashed line) and with Conselice's for- 
mula (dotted line). Even in this case, the Carlberg model 
(m = 3.3) does not match results from simulations, as it 
predicts an increase of the number of mergers with redshift. 



Instead, Conselice's formula is in better agreement with the 
simulated merger rate. In particular, the best-matching pa- 
rameters for equation (0 are 7=5000, m=2.4 and /3=-2.0. 
Then, Conselice's formula reproduces the evolution of the 
CRG formation rate as well as of the total merger rate. 



4 SUMMARY 

We have used cosmological numerical simulations to study 
the rate of mergers of haloes with mass > 5 x 10 12 Mq and 
the CRG formation rate. The large volume combined with 
the selection criterion used to identify halo progenitors al- 
lows us to quantify the cosmic merger fraction, the merger 
rate and the CRG formation rate among haloes. We have 
made comparisons between these and theoretical models and 
the latest available observational data. 

Our main conclusions may be summarized as follows. 

• The merger fraction of progenitors of the present day 
galaxies does not evolve strongly with the redshift between 
Q. 2^7,^.2. Predictions of the merger fraction and merger 
rates are in fair agreement with the current observational 
data, within the great uncertainties. This implies that the 
observed decrease of the cosmic star formation rate since 
z~l is not tied to a disappearing population of major merg- 
ers, at least according to our models. We calculate that 
the number of mergers a progenitor of a halo with mass 
> 5 x 10 12 M Q will undergo from z= 2 to z= 0.2 is N m ~ 2. 
We find that there are still major mergers occurring at red- 
shift lower than z~ 1, mainly due to the late assembly of 
large mass systems. 

• The formation rate of CRGs is a good tracer of the 
merger rate. 

• Assuming that the galaxy interaction rate per unit 
volume is proportional to (1 + z) m , as suggested by pre- 
vious models, we derive m = 3 — 4 from our numerical 
simulations of the concordance cosmological model. How- 
ever, the CRG formation rate as well as the global (ma- 
jor+mmor) merger rate are best-matched by the formula 
iVmcrgor = 7 (1 + z) m exp [,0(1 + z)] (Conselice 2006). 

• The CRG formation rate inferred by simulations is 
in marginal agreement with the observed CRGs between 
0.2 ^z^ 1. However, new redshift measurements are re- 
quired to have a good statistical sample of CRGs. Future 
surveys like zCOSMOS will be able to provide insights on 
the incidence of galaxy mergers and CRGs rates, and will 
be extremely useful to test the rates predicted in the hier- 
archical universe. 
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